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We comprehensively investigate gap solitons and Bloch waves in one-dimensional nonlinear peri- 
odic systems. Our results show that there exists a composition relation between them: Bloch waves 
at either the center or edge of the Brillouin zone are infinite chains composed of fundamental gap 
solitons(FGSs). We argue that such a relation is related to the exact relation between nonlinear 
Bloch waves and nonlinear Wannier functions. With this composition relation, many conclusions 
can be drawn for gap solitons without any computation. For example, for the defocusing nonlinear- 
ity, there are n families of FGS in the nth linear Bloch band gap; for the focusing case, there are 
infinite number of families of FGSs in the semi-infinite gap and other gaps. In addition, the stability 
of gap solitons is analyzed. In literature there are numerical results showing that some FGSs have 
cutoffs on propagation constant (or chemical potential), i.e. these FGSs do not exist for all values 
of propagation constant (or chemical potential) in the linear band gap. We offer an explanation for 
this cutoff. 

PACS numbers: 42.65.-k,42.65.Tg,42.65. Jx, 03.75.Lm 



I. INTRODUCTION 

With the state-of-the-art technology, various nonlinear 
periodic systems have been experimentally realized P, 
131. Typical examples include nonlinear waveguide arrays 
B 01 J optically induced photonic lattices 0] , and Bose- 
Einstein condensates (BECs) in optical lattices [3,0]. In 
these nonlinear periodic systems, there exist two typical 
stationary solutions, Bloch waves and gap solitons, which 
are keys to the understanding of such systems. 

Bloch wave is intrinsic to periodic systems 0- Hence 
the concept of Bloch wave originally introduced for lin- 
ear periodic systems can be straightforwardly extended 
to the nonlinear periodic systems. In both cases, Bloch 
waves are extensive and spread over the whole space Q. 
Nonlincarity, however, will significantly affect the stabil- 
ity of Bloch waves. Whence, the analysis of the stability 
of nonlinear Bloch waves (NBWs) has been a focus of ex- 
tensive research. For example, the nonlinearity-induced 
instabilities of Bloch waves are directly responsible for 
the formation of the train of localized filaments observed 
in various optical systems For another instance, 

the instability of Bloch waves has been experimentally 
observed for BECs in optical lattices [l^, [lj|, where the 
instability is closely related to the breakdown of super- 
fluidity in such systems d, [ill. [T3|. 

In contrast, gap solitons are localized in space and only 
exist in nonlinear periodic systems [J, Q. So far they 
have been found to exist in systems of different natures, 
including nonlinear optical systems P, 0, [20l. [2ll . [2^ . 
BEC systems [ll[ll[2|, even a surface system [2J] . One 
particularly important type of gap solitons is fundamen- 
tal gap solitons (FGSs), whose main peaks are located 
inside a unit cell [11, [H • These FGSs can be viewed 
as building blocks for the higher order gap solitons p^ . 
We note that other types of localized solutions may also 



exist for nonlinear periodic systems, such as gap vortex in 
two-dimensional nonlinear periodic systems [26[. These 
localized solutions persist in discrete models, where they 
are called discrete soliton P, Q and discrete vortex [27j . 
respectively. 

A viewpoint has been floating in the community that 
the NBWs (or higher order gap solitons) can be regarded 
as chains composed of FGSs in one-dimensional nonlin- 
ear periodic systems [H, IH, [2^, [s^l • Such a viewpoint 
would occur to anyone who has observed the almost per- 
fect match between a NBW and the corresponding FGS 
as shown in Fig. [T] In particular, Alexander et al. [2§| 
have found that a new set of stationary solutions, which 
they call gap waves, can be regarded as the intermediate 
states between NBW and FGS. This development is a 
great boost to such viewpoint. However, doubt always 
lingers as people know that a match between a NBW 
and the related FGS can be quite bad for a different set 
of parameters, as shown in Fig. [2lja). Recently, we ap- 
proached this composition relation from a different angle 
[3l|. There, we presumed the existence of this compo- 
sition relation and then investigated how many conclu- 
sions could be drawn from it without any computation. 
These conclusions were eventually verified through nu- 
merical computation. In this way, we were able to make 
a firm claim beyond doubt that the composition relation 
between NBWs and FGSs exists. Note that a similar 
relation was pointed out by Carr et al. for exten- 
sive periodic solutions and solitons for a nonlinear system 
without periodic potential. 

In this paper we explore in detail this composition re- 
lation between NBWs and FGSs We not only offer 
more details on this relationship but also generalize it to 
the systems with focusing nonlincarity. With this rela- 
tion one can draw many conclusions without any com- 
putation. For example, there are n families of FGSs in 
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the nth band gap for the defocusing nonhnearity and the 
FGS of the nth family has n main peaks. All the conclu- 
sions will be discussed in detail and be verified with ex- 
tensive numerical results. Moreover, we have computed 
the Wannier functions from the NBWs [s^] and compared 
them to the FGSs. We find that these nonlinear Wan- 
nier functions match very well with the FGSs. This fact 
seems to suggest that the composition relation between 
NBWs and FGSs is related to the exact relation between 
NBWs and nonlinear Wannier functions. In addition, 
we have analyzed with different methods the stabilities 
of the new-found FGSs and the related gap waves [29j . 
One method is linear stability analysis; the other is the 
so-called nonlinear analysis by integrating the dynamic 
equation with noise (26l | . Our numerical results show that 
not all of these solutions are stable. The stability regions 
are marked out. 

There arc numerical results in literature indicating 
that some FGSs do not exist for all the values of propa- 
gation constant (or chemical potential) in the linear band 
gap [H, [l^, [l3| . Namely, there exists a kind of cutoff. 
Yet it is not clear so far why there is such a cutoff. Here 
we show that there indeed exists such a cutoff arising 
from the mixing of different types of FGSs. This mix- 
ing can be intuitively viewed as a result of a "chemical 
reflection" . 

This paper is organized as follows. In Sec. II, we give a 
brief description of our model equation and show how it is 
related to the concrete systems. In Sec. Ill, we state the 
composition relation between NBWs and FGSs and list 
all the predictions that can be made with this relation. 
We then demonstrate that all the predictions are valid. 
In Sec. IV, the FGSs are compared to nonlinear Wannier 
functions. They resemble each other very well, suggesting 
that the composition relation is related to the well known 
relation between Bloch waves and Wannier functions. In 
Sec. V, the composition relation is applied to construct 
stationary solutions other than Bloch waves with FGSs, 
such as gap waves and multiple periodic solutions. In Sec. 
VI, the stabilities of FGSs and gap waves are examined. 
In Sec. VII, we offer an explanation why FGSs do not 
exist for all values of the propagation constant in the 
linear band gap. Finally, we summarize our results in 
Sec. IX. 



In optics, Eq. ((T|) describes light propagation along the 
z direction in the presence of a periodic modulation in x 
direction. The periodic structure described by V{x)can 
be experimentally realized with waveguide arrays [3| or 
optical inducing technology @ . As routinely used in lit- 
erature, z and X here are respectively scaled to diffractive 
length and beam width. 

In the context of the BEG system, Eq. ([T]) gives the 
description of a BEG in the one-dimensional optical lat- 
tice with z being the time variable. In such case, Eq. ([T|) 
has been scaled as follows: x is in units of A/(27r), z is in 
units of h/{8Erec) and the strength of the optical lattice 
V is in units of SErec with m being the atomic mass, A 
the period of the lattice and Erec ~ ft^7r^/(2mA^) the 
recoil energy. 

For stationary solutions in the form of "^{x^z) = 
(j){x) exp^—ifiz), Eq. ID) is reduced to a z-independcnt 
equation 

--g^ + iycos{x)<j) + a\^\''(b^ ficj). (2) 

Here, fi is referred to as propagation constant in optics. 
Whereas in the BEG system, fi represents the chemical 
potential. In general, NBW and FGS are the two basic 
types of stationary solutions to Eq. For the FGS 

which is localized in space, we can define its norm N as 

|(/>(x)pdx. (3) 

By contrast, the NBW spreads over the whole space, 
hence it is only meaningful to define its averaged norm 
A/" over one period 

(4) 

Jo 

The so defined norm is proportional to the laser strength 
in optics, or the number of atoms for a BEG. 

To avoid confusion, hereafter we will present our re- 
sults and discussions in the framework of optics unless 
otherwise specified. We first consider the defocusing case, 

(7=1. 



II. MODEL EQUATION 

We consider a one-dimensional nonlinear periodic sys- 
tem described by 

^^ = --^ + n-)*+.i*pv,, (1) 

with V{x) being a periodic function. Without loss of 
generality, we will use V{x) = v cos{x) throughout this 
paper. The a in Eq. ([T]) indicates the type of nonhn- 
earity: (7 = 1 for the defocusing (or repulsive) case and 
cr = — 1 for the focusing (attractive) case. 



III. COMPOSITION RELATION IN THE 
DEFOCUSING CASE 

Without the nonhnear term, Eq. ^ is the well- 
known Mathieu equation (ssj . Its physical solutions are 
Bloch waves defined hy 6n,k(x) = exp(ikx)il'n,kix) with 
4'n,k{x) = ipn,k{x + 2'K) Here k is the Bloch wave vec- 
tor and n is the band index; the /x„(fc) form Bloch bands 
as k varies through the Brillouin zone (BZ). There exist 
band gaps between different Bloch bands indexed by n, 
where the physical solutions are forbidden. 

With the addition of the nonlinear term, the physical 
solutions of Eq. ^ become admissible in the linear band 
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gaps. One such typical solution is the gap soliton. Since 
the propagation constants of gap soliton only take values 
inside the linear band gaps, no linear counterpart exists 
for gap soliton. Among various gap solitons, there is 
a particularly important class called FGSs, whose main 
peaks locate inside one unit cell. 
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FIG. 1: NBWs (dotted lines) of the first nonlinear band and 
FGSs (solid lines) in the first linear band gap for v = f .5 and 
^ = -0.3. NBWin(a) is at the center of BZ with A/" = f.6908; 
NBW in (b) is at the edge of the BZ with N = 1.6738. 




FIG. 2: NBWs (dotted lines) of the first nonlinear band and 
FGSs (solid lines) in the first linear band gap for v = 1.5. (a) 
NBW is at the center of BZ with M = 0.0027 and p. = -0.92 
is near the top of the first linear Bloch band; (b) NBW is at 
the edge of the BZ with M = 3.2194 and [i = 0.151 is close to 
the bottom of the second linear Bloch band. 

Despite the nonlinearity, Eq. ^ still admits the Bloch 
wave solutions. However, the nonlinear term will modify 
Bloch waves and Bloch bands. For the defocusing case, 
the nonlinearity not only changes the shapes of the Bloch 
bands but also moves them up into the linear band gap 
[sl ]. Here, the strength of nonlinearity is completely de- 
scribed by the norm M . If M is lowered to zero, these 
nonlinear bands will move down and be reduced to the 
linear Bloch bands. When M is increased, bands will 
move up continuously without limit. 

It is clear from the above discussions that there exist 
two types of solutions to Eq. 0, NBW and FGS, for a 



given [L in the linear band gap. In Fig. [1] we have plot- 
ted both NBWs and FGSs for [i = -0.3 in the middle 
of the first linear band gap. Fig. HUa) (or (b)) is for a 
NBW at the center (or edge) of the BZ. In this figure, 
a nearly perfect match is found between the NBWs and 
corresponding FGSs inside one unit cell. These numer- 
ical results therefore strongly suggest that FGS can be 
considered as the building blocks for NBWs at either the 
center or the edge of the BZ. In other words, NBW at 
the center of the BZ can be viewed as an infinite chain 
composed of FGSs while NBW at the edge of BZ is built 
by FGSs with alternative signs. 

However, such perfect match does not exist for all \i 
in the linear band gaps. After checking various values of 
/.t, we find that the match between the NBW and FGS 
is very good except in a narrow region near the edge of 
the linear bands. Two typical results are shown in Fig. 
m Fig- IHa) is for the case of [l near the top of the first 
Bloch band while Fig. [l{b) is for the case of [i close 
to the bottom of the second Bloch band. It is evident 
from Fig. [2ja) that for a [i close to the edge of the first 
linear band, the NBW and the FGS does not match well. 
This mismatch casts strong doubt on the validity of the 
claim that a NBW can be regarded as an infinite chain 
composed of FGSs. 

In this work, the gap solitons are numerically obtained 
by using the relaxation method in the coordinate space 
(itI [ist while the Bloch waves are numerically foun d by 
applying the relaxation method in the Fourier space [3J] . 



Direct predictions from the composition 
relation 



We now take a different view at the above observed 
composition relation between the FGSs and NBWs. We 
shall first presume the existence of this composition re- 
lation and then try to draw as many conclusions as pos- 
sible. By verifying these conclusions, we will justify this 
composition relation a 'posteriori. In this sprit, the fol- 
lowing predictions can be immediately drawn without 
any computation. 

(1) For defocusing nonlinearity, there is no FGS in the 
semi-infinite linear gap below the lowest Bloch band. 

As defocusing nonlinearity can be regarded as a re- 
sult of repulsive interaction, its addition to the system 
will increase the system energy and therefore move the 
nonlinear Bloch bands up relative to their linear coun- 
terparts. This means that there is no NBW for the /i 
in the semi-infinite gap. According to the composition 
relation, one can then conclude that there is no FGS in 
the semi-infinite band gap. 

(2) There exist n different families of FGSs in the nth 
linear band gap for defocusing nonlinearity. 

In order to show this, we have plotted the linear and 
nonlinear Bloch bands in Fig. [31 Note that all the NBWs 
in the same nonlinear Bloch band share the same nonlin- 
earity Af. As already discussed above, the nonlinearity Af 
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FIG. 3: Linear and nonlinear (defocusing case) Bloch bands 
for V = 1.5. Solid lines are linear bands, LBi represents the 
ith linear band. The widthes of linear bands are marked by 
shadow areas, (a) The first nonlinear Bloch band (dotted 
line) in the first linear band gap with TV = 1.6908; (b) the 
first (dotted) and second (dashed line) nonlinear bands in 
the second linear band gap with TV = 4.8437 and 1.5268 re- 
spectively; (c) the first (dotted), second (dashed), and third 
(dash-dotted line) nonlinear bands in the third linear band 
gap with TV = 8.3478, 4.7363, and 1.0242 respectively. 



can move the nonlinear Bloch bands up. Therefore, the 
mth nonlinear Bloch band can be lifted into the n > m 
linear band gaps. For example, as Af increases, the first 
nonlinear Bloch band can be lifted into the first, second, 
third, and all other linear band gaps. This implies that 
there exists only one nonlinear Bloch band ( the first non- 
linear Bloch band) in the first linear band gap [FiglSja)]; 
two nonlinear Bloch bands (the first and second) in the 
second linear band gap [FiglSlJb)]; three nonlinear Bloch 
bands (the first, second, and third) in the third linear 
band gap [Fig[3ljc)]; and so on. In other words, there are 
n different NBWs in the nth linear band gap. NBWs in 
the different Bloch bands have different characters. With 
the composition relation, one can immediately conclude 
that there are n different families of FGSs in the nth 
linear band gap. 

In Ref. this rising nonlinear Bloch band by non- 
linearity was noted to be useful for analyzing gap waves 
(or truncated Bloch waves). 

(3) In the nth linear band gap, the mth (m < n) family 
of FGSs exists only above a threshold value of norm N 
whereas the nth family does not have such a value. 

Generally, one must increase nonlinearity Af over a 
critical value to move the mth nonlinear Bloch band up 
into the nth (n > to) linear band gap while there is no 
such a critical value to lift the nth nonlinear Bloch band 
into the nth linear band gap. An example is shown in 
Fig. injb) . In order to lift the first nonlinear band into the 
second linear band gap, nonlinearity TV must be beyond 
a threshold value while there is no such a value to move 
the second nonlinear band into the second linear band 
gap. This analysis, combined with the composition rela- 
tion, leads us to predict that there is a threshold value 
of norm TV for the mth {m < n) family of FGSs in the 
nth linear band gap while the nth family have no such a 
threshold value. 



(4) The nth family of FGSs has n main peaks inside 
one unit cell ( or an individual well in the periodic poten- 
tial). 

The linear Bloch waves in the nth linear Bloch band 
originate from the nth bound state of an individual well 
of periodic potential. Since the nth bound state has n — 1 
nodes, the linear Bloch waves have n main peaks in one 
unit cell. This character is shared by the Bloch waves 
belonging to the nth nonlinear Bloch band. Therefore, 
as the building blocks of NBWs, the nth family of FGSs 
should have n main peaks. 
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FIG. 4: NBWs (dotted lines) and FGSs (solid lines) in the 
second linear band gap for u — 1.5 and fi — 0.6. (a) A first- 
family FGS and its corresponding NBW of the first Bloch 
band at the center of the BZ with TV = 5.0598; (b) a first- 
family FGS and its corresponding NBW of the first Bloch 
band at the edge of the BZ with TV = 4.8443; (c) a second- 
family FGS and its corresponding NBW of the second Bloch 
band at the center of the BZ with TV — 1.5268; (d) a second- 
family FGS and its corresponding NBW of the second Bloch 
band at the edge of the BZ with TV = 1.7417. 



B. Verification of the predictions 

In the following, we check the validity of the predic- 
tions listed above. 

(1) The first prediction is consistent with the well- 
known and extensively proved fact that gap solitons do 
not exist in the semi-infinite gap for defocusing nonlinear- 
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FIG. 5: NBWs (dotted lines) and FGSs (solid line) in the 
third linear band gap for v = 1.5 and fi = 1.4. The solid 
lines are for FGSs and dotted lines for NBWs. (a,b) A first- 
family FGS; (c,d) a second-family FGS; (e,f) a third-family 
FGS. (a,c,e) NBWs at the center of the BZ in the first, sec- 
ond, and third bands with Af = 9.2162, 4.7363, and 2.052, 
respectively; (b, d, f) NBWs at the edge of the BZ in the 
first, second, and third bands with Af = 8.3478, 5.6197, and 
1.0242, respectively. 



ity [l5|,ll6|. As a result, this prediction can be considered 
as the confirmation of a known result. 

(2) We resort to the numerical computation to verify 
the second prediction. As it is impossible to exam every 
linear band gaps, we foucus on the second and third linear 
band gaps. We indeed find two families of FGSs in the 
second linear band gap and three families of FGSs in the 
third linear band gap. They are shown and compared to 
the corresponding NBWs in Figs. IHandO 




FIG. 6: (Color onhne) Norms of FGSs and NBWs as func- 
tion of the for v = 1.5. Shaded areas are linear bands. 
Dotted(red), dashed(green), solid(black) lines represent Bloch 
waves at the BZ center, Bloch waves at the BZ edge and FGSs 
respectively. 



We note that the second family of FGSs are called 
subfundamental gap solitons in literature This in- 

dicates that people were not expecting other FGSs to be 
found. In other words, the existence of the third family 
of FGSs as shown in Figs. [5] (e) and (f) is a surprise to 
many. Our results here also show that all the families 
of FGSs should be regarded equally fundamental as the 
corresponding NBWs are equally important. 

(3) In Fig. [21 we have plotted the TV as a function 
of the /i for three different families of FGSs. It is clear 
from the figure that in the second linear band gap the 
first family of FGSs exists only when their N > 3.4140 
while the second family exists for arbitrary small norm, 
having no threshold value. In the third linear band gap, 
the first and second family exist only when their norms 
N > 8.6105 and N > 5.0152, respectively. In contrast, 
the third family of FGSs has no such threshold value. 
For comparison, we have also plotted the norms M of the 
corresponding NBWs versus fi in Fig. [6l These NBW 
norms N' match quite well with the FGS norms A''. 

(4) We have also found that the number of main peaks 
of a FGS in a well is just what we have expected. In 
order to demonstrate this clearly, we have plotted Fig. [7| 
the three families of FGSs in the third linear band gap 
along with the periodic potential. It is clear that the first 
family of FGSs has one main peak in an individual well 
[Fig. EJa)]; the second and third families of FGSs have 
two and three main peaks in a unit cell, respectively, as 
shown in Figs. [Tl^b) and (c). 



IV. NONLINEAR WANNIER FUNCTION AND 
FUNDAMENTAL GAP SOLITON 

For a linear periodic system, there is another well- 
known function, Wannier function, which is localized in 
space. In one dimension, the Wannier function Wn{x) is 
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where Xj is the location of the j th well in the periodic po- 
tential and the summation runs over all potential wells. 
The relation of Eq. ^ is still valid for nonlinear pe- 
riodic systems and the corresponding Wannier function 
can be called nonlinear Wannier function [s^ . At either 
the center or the edge of the BZ, Eq. ^ becomes 



(6) 



FIG. 7: Typical profiles of FGSs in the third linear band gap 
for v = 1.5, and fi = 1.4. (a) The first family; (b) the second 
family; (c) the third family. Dotted lines in the bottom of 
each figure mimic the periodic potential. 
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where + is for the center and — is for the edge. This re- 
lation is very similar to the composition relation between 
NBWs and FGSs that we have just established. Hence, 
a question naturally arises whether the FGSs bear any 
relation to the Wannier functions. 

In Fig. [5l wc have plotted nonlinear Wannier functions 
and the corresponding FGSs for the first Bloch band to- 
gether for comparison. The Wannier functions are com- 
puted from the Bloch waves in a standard way [3^ . For 
the nonlinearities M considered here, the first nonlinear 
Bloch band lies completely in the first linear band gap. 
As a result, there are two different FGSs for this band, 
one corresponding to the NBW at the BZ center and the 
other to the NBW at the BZ edge. Both of the FGSs are 
plotted in Fig. [5] and are found to match the Wannier 
functions very well. In fact, the match gets better as the 
periodic potential gets stronger. Since the Wannier func- 
tion is normalized, we have scaled them by a factor ^pN 
for comparison with the FGSs in Fig. [5l 

The excellent match between the FGSs and the Wan- 
nier functions suggests that they are related. This rela- 
tionship may be intuitively understood in the following 
way. Although a Bloch wave is a solution of a periodic 
system, the Wannier function is not. In a linear periodic 
system, any localized wave function, including a Wannier 
function, will spread in space. In a nonlinear periodic sys- 
tem, it seems that the Wannier function can be modified 
slightly and become a solution of the system in the form 
of a FGS. 



V. FUNDAMENTAL GAP SOLITONS: 
BUILDING BLOCKS FOR STATIONARY 
SOLUTIONS 



FIG. 8; (Color online) Nonlinear Wannier functions for the 
first nonlinear Bloch band and the corresponding FGSs. Solid 
(black) lines are for Wannier functions, dashed (green) lines 
are for the FGSs corresponding to the NBWs at the band 
edge, and dotted (red) lines are for the FGSs corresponding 
to the NBWS at the band center, (a) v = 0.2 and TV = 0.5027; 
(b) V = 0.4 and N = 0.7540; (c) v = 1.5 and N = 1.6908. 



The composition relation between NBWs and FGSs 
suggests that the FGSs are really fundamental and can be 
viewed as building blocks for other stationary solutions 
to a nonlinear periodic system, such as high-order gap 
solitons. Our numerical results and existing results in 
literature [H, [H, [H, [s^l fully support this view. In the 
following, we show a few examples. 



related to the Bloch wave c/)„.fe as follows Q 



A. Gap waves 

In the first example, we construct new solutions by 
putting finite number of FGSs together. Such solutions 
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FIG. 9: Gap waves (solid lines) and the corresponding NBWs 
(dotted lines) in the third linear gap for f — 1.5 and fi = 1.4. 
The NBWs in (a), (c), and (e) are at the BZ center with 
Af = 9.2162, 4.7363, and 2.052 respectively; the NBWs at the 
BZ edge in (b), (d), and (f) have JV = 8.3478,5.6197, and 
1.0242 respectively. 



are called high order gap solitons in Rcfs. [l8|,|30l. There 
are numerous ways to build these high order gap solitons. 
For instance, one can use FGSs from different families 
and put them together with either the same phase or the 
opposite phase. Here we mainly focus on a class of high 
order gap solitons called gap waves by Alexander et al. 
since they can be viewed as truncated NBWs [1^. For 
gap waves, all the constituent FGSs come from the same 
family. These gap waves can be viewed as intermediate 
states between NBWs and FGSs. With one FGSs, two 
different classes of gap waves can be constructed. The 
first class of gap waves, which we call GW-I, are built 
by putting FGSs side by side similar to the NBW at the 
center of the BZ. The second class, called GW-II, are 
composed of FGSs pieced together with alternative signs 
similar to the NBW at the BZ edge. Some typical gap 
waves in the third linear band gap and the corresponding 
NBWs are plotted in Fig. M Figs. M,i^), (c), and (e) are 
GW-I composed of the first, second, and third families 
of FGSs respectively; Figs. Mh), (d), and (f) are GW-II. 
Note that gap waves and the corresponding NBWs have 
the same propagation constant ^. 
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FIG. 10: Gap waves composed of different number of FGSs 
for !^ = 1.5 and fi = —0.92. Solid lines are for gap waves; 
dotted hne is for the NBW with TV = 0.0027. 



We have also plotted gap waves composed of different 
number of FGSs in Fig. \W\ along with the corresponding 
NBW. The parameters in this figure are the same with 
ones in Fig. ^a), where the mismatch between the FGS 
and the NBW is obvious. Fig. [10] shows an interesting 
trend that the match between the gap waves and the 
NBW improves as the number of FGSs increases. This 
can be viewed as another supporting evidence for the 
composition relation. 



B. Multiple periodic solutions 

Multiple periodic solutions are extensive states like 
Bloch waves, but with multiple periods. Mathemati- 
cally, they are defined as (j){x) ~ cxx)(ikx)'ipk(x) with 
ipk{x) = ipkix + 2p7r) and p > 1 being an integer [3q . 

The composition relation between FGSs and NBWs 
can be generalized to construct multiple periodic solu- 
tions. Here we show two typical kinds of multiple peri- 
odic solutions in Fig. [TTJ one is double periodic solution 
[Fig. IUJa)]; the other is triple periodic solution [Fig. 
[TlT b)]. The double periodic solution is constructed with 

the pattern | " and the triple periodic solution is 

built with the pattern "++H " . One can build other 

multiple periodic solutions with other patterns. In fact, 
one can also build the multiple periodic solutions with 
FGSs from different families. The odd-periodic solutions 
were speculated to exist in Ref. [36| . 



C. Composition relation in the presence of a loop 

We have seen in Fig. [3] that all the nonlinear bands 
move up with the increasing defocusing nonlinearity. It 
is known that a more dramatic change can happen when 
the nonlinearity is large enough: loop structures emerge 
at the BZ edge for the first band and at the BZ center 




FIG. 11: Multiple periodic solutions for v = 1.5, \i = —0.1. 
(a) period doubled solution; (b) triple periodic solution. Dot- 
ted lines are multiple periodic solutions, solid lines are cor- 
responding FGS. Lines in the bottom of each figure mimic 
periodic potential. 



for the second band [3J, |37[. The critical value of the 
nonlinearity for the loop to appear \b M > ^-kv [s^I . Our 
studies show that the composition relation still holds in 
the presence of such a loop. As an example, we demon- 
strate this relation for a NBW sitting at the loop top 
of the first band in Fig. [121 In this figure, the looped 
nonlinear Bloch band is already in the third linear band 
gap. 



VI. STABILITY PROPERTIES OF 
FUNDAMENTAL GAP SOLITONS AND GAP 
WAVES 

A question that must be asked for a stationary solu- 
tion (or a fixed point solution) in a nonlinear system is 
whether the solution is stable. The unstable solution is 
very sensitive to small perturbations. The stabilities of 
the NBWs have been discussed thoroughly in Rcf. 8, 3^ . 
In the following, we shall focus only on FGSs and gap 
waves. 

We use two different approaches to examine the stabil- 
ities of the stationary solutions. The first is called linear 
stability analysis. It is done by adding a perturbation 
term to a known solution 



5'(z,a;)= -f A0(z, x) exp(— i/iz), 



(7) 




1.00 



120 



FIG. 12: (a) Loop structure of the first nonlinear band lying 
in the third linear band gap for v = \ and J\f /2n = 1.1344. 
Shadow areas are the linear bands and the blank area is the 
third band gap; (b) NBW (dotted line) for the solid point 
in (a) with = 1.3 and the corresponding gap wave (dashed 
line) and FGS (solid line) in the third band gap. 



perturbation and 4){x) is the stationary solution. Plug- 
ging Eq. ([7]) into Eq. ([1]) and keeping only the linear 
terms, we obtain as follows 
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where A(/)(z, x) — u{x) cxp{i6z)+w*{x) CKp{—iS*z) is the 



FIG. 13: Stability of the second family FGS in the second 
linear gap. (a) Maximum imaginary part of 5 as a function of 
fi; (b) the eigenvalue plane of 5 at /i = 0.6. The inset is the 
profile of the corresponding gap soliton. i/ = 1.5. 
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with £ = - iycos{x) - 2|0p + /i. In Eq. if the 

eigenvalue S has imaginary parts, the solution of (t>{x) is 
unstable; otherwise, the solution is stable. 

In the second method, the perturbed solution in Eq. 
([7]) is used as the initial condition for Eq. (jT]). Its evolu- 
tion is then monitored numerically. If its deviation from 
grows as the system evolves, the solution (j){x) is 
unstable; it is stable otherwise. The stability checked by 
this method is called nonlinear stability |26| . 

A. Fundamental gap solitons 



The above results show that the newly-found third 
family FGS can be stable and therefore should be ob- 
servable in experiment. Usually in order to observed gap 
solitons experimentally, initial input beam profile should 
be close to the desired soliton profiles [H, H^. We pro- 
pose to observe the third family FGS using two local- 
ized laser beams, whose wavelength is much shorter than 
the period of a waveguide, to form an interference pat- 
tern with three large peaks in a unit cell of the periodic 
waveguide. 



Our linear stability analysis shows that the first family 
of FGSs in the first and second band gaps are stable con- 
sistent with Ref. [iBl- However, they become unstable in 
a small area near the band edges in the third gap. 
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FIG. 14: Stability of the third family of FGSs in the third 
gap. (a) Maximum imaginary part of 5; (b) the eigenvalue 
plane of 5 at /i = 1.23. The inset is the corresponding gap 
soliton. 1/ = 1.5. 

The stability of the second family of FGSs in the sec- 
ond linear band gap is shown in Fig. Iisr a). where the 
stability is measured by the maximum value of the imag- 
inary part of the eigenvalues S. If the maximum value is 
zero, the solution is stable; otherwise, it is unstable. It is 
clear from Fig. [TSlJa) that this family of FGSs is stable 
when their fi are smaller than a critical value near 0.5. 
For other values of /i above this critical value, the soli- 
tons become unstable. Fig. [TSTb) are the eigenvalues for 
a soliton illustrated in the inset. We see that the eigen- 
values 5 are mostly real and become complex only in a 
very small region. The stability of the third family in 
the third linear band gap is shown in Fig. [M) The result 
is very similar to that of Fig. [13] except that the stable 
region is much smaller. Fig. I14r b) shows an example of 
the eigenvalue plane. 

To double check the stability results, we have prop- 
agated perturbed FGSs by numerically solving Eq. ^ 
with the split-step Fourier method. Gaussian distributed 
random noise is added to FGSs for the initial wave func- 
tion. The propagation results shown in Fig. [15] agree 
with our linear stability analysis. Fig. fTSTa) is the prop- 
agation of a stable solution while Fig. ITSl b) is for an 
unstable solution. 
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FIG. 15: Propagations of gap solitons in the presence of Gaus- 
sian distributed random noise with variance ct'^ = 0.01. (a) 
Evolution of a third family FGS at fi = 1.21 which is sta- 
ble as indicated by the linear analysis results in Fig ll4r a): 
(b) evolution of an unstable third family FGS at /i = 1.23, 
which is unstable as indicated by the linear analysis results 
in FigdHa) (Max[Imag((5)]=0.0115). 



B. Gap waves 

The stability of gap waves is also analyzed. In general, 
the GW-I's corresponding to the first band are stable 
in the first and second gaps but unstable in the third 
gap. The stability of gap waves composed of the second 
family of FGSs in the second gap is shown in Fig. [TH] 
These gap waves contain either 4, 6, or 9 FGSs. These 
gap waves are characterized by the norms in Fig. fTHT a). 
Fig. flGT b) demonstrates that GW-I are always unstable. 
GW-II are stable in a small regime near the top of the 
second band, but they are unstable for other values of 
/i as shown in Fig. [TU^c). The propagation of the gap 
waves with noise confirms our stability analysis as shown 
in Fig. [T7] As the three curves fall almost on top of 
each other in Figs. riHT b.c). we find that the stability of 
gap waves are independent of how many FGSs they have. 
Our analysis shows that gap waves other than the types 
discussed above are unstable. 
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FIG. 16: (Color online) Stabilities of gap waves composed of 
the second family of FGSs in the second gap. (a) the families 
of gap waves. Circles are for GW-I while triangles for GW-II. 
Shadow areas are the linear bands; (b) maximum imaginary 
part of S for GW-I; (c) maximum imaginary part of 5 for 
GW-II. 




FIG. 18: Norms of fundamental gap solitons as functions of fi 
for focusing nonlinearity. Shaded areas are for linear bands; 
lank areas are for gaps, f = 1.5. 




FIG. 17: Propagations of gap waves in the presence of a 
Gaussian distributed random noise with variance = 0.01. 
(a) propagation of GW-II containing six FGSs ai jj. — 0.24, 
which is stable indicated by the linear analysis result shown 
in Fig llBf cl ; (b) propagation of GW-II with six second family 
FGSs at /I = 0.26, which is unstable as shown in Fig |16l (c) 
(Max[Imag((5)]=0.0162). Note the different scale of z in (a) 
and (b). 



VII. COMPOSITION RELATION IN THE 
FOCUSING CASE 

We have concentrated on the defocusing case. We 
now turn to the focusing case. For focusing nonhnearity, 
Bloch bands and NBWs still exist. Unlike the defocusing 
case, focusing nonlinearity causes nonlinear bands move 
down. As a result, the predictions made from the com- 
position relation are different from the defocusing case. 
(1) There exist infinite number of families of FGSs in the 
semi- infinite and finite linear band gaps. It is because 
infinite number of bands can move into a given linear 
band gap with increasing focusing nonlinearity. (2) In 



the nth linear gap {n > with n, ~ for the semi-infinite 
gap), the nth family and other lower order fmailies of 
FGSs do not exist. (3) In the nth linear gap, only the 
(n + l)th family FGSs exist for an arbitrary small values 
of norm while all other families of FGSs exist only for 
norms above certain threshold values. 

These predictions are confirmed by our numerical com- 
putation for the first three bands and the corresponding 
three band gaps. The results are summarized in FigllSl 
where the norms of different FGSs are plotted as func- 
tions of fi. As shown in this figure, corresponding to these 
three bands, there are three families of FGSs in the semi- 
infinite linear band gap, two families of FGSs in the first 
linear band gap, and one family of FGSs in the second 
band gap. In other words, there exist no first family of 
FGSs in the first linear gap and there exist no first and 
second families of FGSs in the second gap. Another fea- 
ture in the figure is the threshold values of norm for some 
families of FGSs. In the semi-infinite gap, the threshold 
value for the second family of FGSs is 3.2824 and for the 
third family 7.2940. In contrast, the first family have 
no such threshold value. In the first gap, the threshold 
value of norm for the third family of FGSs is 4.0616 while 
the second family has no threshold value. The match be- 
tween the norms of NBWs and FGSs is similar to that in 
the defocusing case. 

Combining with the results for the defocusing case, we 
have an interesting observation: in the nth linear Bloch 
band gap, the first n families of FGSs exist for the de- 
focusing nonlinearity and the other families ((n -I- l)th, 
(n + 2)th, • • • ) exist for the focusing nonlinearity. 



VIII. CHEMICAL REFLECTION OF FGS 

One may have noticed in Figs. [H] and [T5] that some 
families of FGSs do not exist for all values of the /i in the 
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FIG. 19: (Color online) Chemical reflection of FGSs into high 
order gap solitons in the (A'', /x) plane for defocusing case with 
v = 1.5. (a) The N-fj, curves for the first family of FGSs 
in the second gap (solid-red line) and the higher order soli- 
tons (dashed-green line) composed of one first-family FGS 
and two second-family FGSs; (b) the A'^-^ curves for the first 
family of FGSs in the third gap (solid-red line) and the soli- 
tons (dashed-green line) consisting of one first-family FGS 
and two third-family FGSs; (c) the A'^-^ curves for the second 
family of FGSs (solid-red line) and the solitons (dashed-green 
line) composed of one second-family FGS and two third-family 
FGSs. (d, e, f) soliton profiles corresponding to the labeled 
points in (a, b, c), respectively. Green lines are for profiles 
of high order solitons, red lines are FGSs, and black lines are 
profiles of the second family FGS in (d) and the third family 
FGS in (e) and (f). The lines in the bottom represent periodic 
potential. 



linear band gaps. For example, in Fig. [6l the second and 
third families of FGSs in the third linear band gap do 
not exist for n near the edge of third linear Bloch band. 
The N-fj, curves for these two families end at ^ = 1.2658 
and /i = 1.2612, respectively, which are away from the 
right edge of third linear band at = 1.1422. This cut- 
off phenomenon was noticed before [1^. However, to our 
best knowledge, no one is sure why this cut-off happens. 
In the following, we show that this cut-off is caused by 
the mixing of different types of FGSs, which can be in- 
tuitively viewed as a result of a "chemical reflection" . It 
will be discussed for both defocusing and focusing non- 
linearities. 

We consider first the defocusing case. We have re- 
plotted the N — ^ curves (solid lines) in Figs. fToya.b.c). 
where the cut-offs exist. In these three figures, we have 
also plotted the N-fi curves (dashed lines) for three differ- 
ent classes of high order gap solitons. Interestingly, they 
are connected smoothly to the curves for the FGSs. As 
illustrated in Figs fTW d.e.f). we find after careful exam- 
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FIG. 20: Chemical refiection of the second and third families 
of FGSs with higher order gap soliton in the semi-infinite gap 
for focusing case, f — 1.5. (a) the third family of FGSs (solid- 
red line) and higher order gap solitons (dashed-green line) 
composted of one third- family FGS and two first-family FGSs. 
(b) the second family of FGSs (solid-red line) and higher order 
gap solitons (dashed-green line) composted by one second- 
family FGS and two first-family FGSs. Dotted lines in (a) 
and (b) are the first linear band, (c) and (d) soliton profiles 
corresponding to labeled point in (a) and (b) respectively. 
Green lines are for profiles of high order solitons, red lines are 
FGSs, and black lines are profiles of the first family FGS. The 
lines in the bottom represent periodic potential. 



ination that the high order gap solitons for the dashed 
curve in FigfTWa) are composed of a first-family FGS 
sitting in one site and two second-family FGSs sitting in 
two neighboring sites, the high order solitons in Fig fTOT b) 
composed of a first-family FGS and two third-family 
FGSs, and the high order solitons in FigfTWc) consists 
of a second- family FGS and two third-family FGSs. 

To help us understand the turning N — fi curves, we 
have developed an intuitive picture to visualize this re- 
sult. We use the case in Fig. [191 a) as an example. If one 
imagines an "atom" moving along the lower curve for 
the FGSs in the (-/V, /x) plane in Fig. [TWa). this "atom" 
gets reflected by the "repulsive walls" of the second lin- 
ear band. Moreover, a "chemical reaction" occurs during 
the collision between this "atom" of FGS and the "wall" , 
which may be viewed as a crystal made of "atoms" of the 
second- family FGS. The result of this reaction is that the 
"atom" changes its nature to a "molecule" of high order 
soliton by picking up two second-family FGSs from the 
second linear band. This "chemical reaction" similarly 
occurs in Fig. llQr b.c). Based on this intuitive picture, 
we call this cut-off phenomenon in Figs. fTQla.b.c) chem- 
ical reflection. 

Note that a similar turning N — fi curve was also found 
for gap vortexes in Ref. [4l| and gap waves in Ref. [s^ . 

The cut-off phenomenon for the focusing case as shown 
in Fig. [THlcan be similarly be viewed as the result of the 
chemical reflection. In the focusing case, as shown in 
Fig. [THl except the lowest family of FGSs in each linear 
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gap, all other families have cut-offs in the propagation 
constant The cut-off phenomenon can be similarly 
viewed as the result of the chemical reflection as demon- 
strated in Fig. [201 Fig. is for the third family of 
FGSs in the semi-infinite gap, whose N-jj curve is found 
to be connected smoothly to a class of high order gap 
solitons composed of a third family FGS in one site with 
two FGSs of the first family in its neighboring sites [see 
Fig. [20t c')] . The case for the second family of FGSs in the 
semi-infinite gap is shown in Fig. I20f b). where the high 
order solitons consist of one FGS of the second family 
and two FGSs of the first family [see Fig. fIU\d)]. 

IX. CONCLUSION 

In conclusion, we have demonstrated that a composi- 
tion relation exists between FGSs and NBWs for both 



defocusing and focusing nonlinearities. Based on the 
composition relation, we have drawn many conclusions 
about the properties of FGSs directly from Bloch band 
gap structures without any computation. All the predic- 
tions have been examined and confirmed. All our studies 
point to one important conclusion that the FGSs are re- 
ally fundamental and they serve as building blocks for 
other stationary solutions in one-dimensional nonlinear 
periodic systems. 
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